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Within the framework of exact quantum electrodynamics in dielectric, we study the topological 
Berry phase of a classically pumped A-type three-level atom, prepared initially in a superposition 
of its two pumped levels and located near a planar dielectric half-space with model permittivity 
of Lorentz type, interacs with a vacuum field. The effects of material losses, expressed in terms 
of the Green tensor of the dielectric-matter formation including dispersion and absorption, on the 
topological Berry phase of the photon transition between the pumped levels have been studied. The 
outcomes are compared with that of an atom in free space. We expect, the additional noise due to 
presence of the dielectric will modify, in a unified way, the behavior and values measured of the phase 
shift between the wavefunction evolution and its original state. It is shown that for small separation 
distance between atom and the dielectric surface, the angle separation between the wavefunction 
evolution and the original state reduced noticeably that is of strong applications in construction of 
the universal quantum logic gates. Further, the study opens routes for new challenging applications 
of atomically systems as various sources of coherent light emitted by pumped atoms in dielectric 
surroundings. 
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I. INTRODUCTION 

In experiments, the use of optical instruments, which are macroscopic material bodies, such as beam splitters or 
cavities, requires careful examination with regard to their action on the light under study. In principle, such bodies 
could be included as a part of the matter to which the radiation field is coupled and treated microscopically [lj. 
Dielectric matter plays an important role in optics, because (passive) optical instruments are typically composed 
of dielectrics. In quantum optics an important consideration is the influence of the presence of material bodies on 
the quantum statistics of the light. Commonly, dielectric matter is characterized by the permittivity; a complex 
function of frequency with real part (responsible for dispersion) and imaginary part (responsible for absorption), 
which describes the response of the matter to the electric field. From statistical mechanics it is clear that dissipation 
is unavoidably connected with the appearance of a random force which gives rise to an additional noise source of 
the electromagnetic field. Hence, any quantum theory that is based on the assumption of a real permittivity can 
only be valid for narrow-bandwidth fields far from medium resonances where absorption can safely be disregarded. 
In 1946, Purcell [2|, through Purcell effect, has stimulated series of theoretical and experimental works . The effect 
lies in the fact that the spontaneous decay rate of an excited atom essentially depends on whether or not the atom 
is located near optical inhomogeneities and interfaces of media with differing optical properties. In fact, spontaneous 
decay, which is a property of spontaneous emission of an excited atom, depends sensitively on the photonic spectral 
density of states that are involved in the atomic transition at a chosen location of the atom. However, depending 
on the specific configuration of inhomogeneities (interfaces) ; expressed in terms of the Green tensor of the dielectric- 
matter formation including absorption and dispersion, the atomic spontaneous decay rate may both increase and 
decrease compared with that of the same atom in free space according to the modifications in the photonic density 
of states due to the presence of macroscopic bodies. It is well known that the spontaneous decay of an excited atom 
can be strongly modified when it is placed inside a microcavity 0, 0j- Purcell effect took on special significance 
in view of rapid progress in physics of low-dimensional nanostructures. It was shown to be of great importance 
for microcavities Q, optical fibers @, photonic crystals Q, semiconductor quantum dots [8]. At this point, a very 
interesting question is that of the topological phase (i.e., the adiabatic geometric phase) of the wavefunction due to 
the presence of the matter. Obviousely, the quantum vacuum effects, such as, the Casimir effect (HQ, anomalous 
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magnetic moment of electron, vacuum polarization leading to Lambs shift as well as Casimir-Polder potentials 
[J EI) atomic spontaneous radiation and its modification in photonic crystals 12], magnetoelectric birefringences 
of quantum vacuum fl3l | and vacuum contribution to the momentum of anisotropic media (e.g., magnetoelectric 
materials) 14], will strongly influenced due to the presence of these dielectric bodies and as a consequence the change 
in phase shift between the evolution of the wavefunction and its original state. 

Geometric phases have received special attention of many researchers since Berry showed that there exists a topo- 
logical phase in quantum-mechanical wavefunction in adiabatic quantum process where the cyclic evolution of wave- 
function yields the original state plus a phase shift, which is a sum of a dynamical phase and a geometric phase shift 
Berry phase [lj| or geometric phase, which does not have classical correspondence, became a focus point in 
modern physics. It describes a phase factor gained by the wavefunction after the system undergoes an adiabatic and 
cyclic evolution, which reflects the topological properties [iH [I?} of the state space of the system and has untrivial 
connections with the character of the system [l8[ , especially with the entanglement pjl, [2(| • The associated topolog- 



ical Berrys phase is known as the Pancharatnam phase |21j , who introduced the concept of a geometric phase during 



his studies of interference effects of polarized light waves. Recently, the Berry phase was introduced into quantum 
computation to construct a universal quantum logic gates that may be robust to certain kinds of errors dm}, it was 
shown [23|, [27l - l36j that the geometric phase shift can be used for generation of built-in fault-tolerance phase shift gates 
in quantum computation. Chen et al., [371 ] investigated the geometric phase gate in a system with many identical 
three-level atoms confined in a cavity under decay, driven by a classical field. They (37| proposed an efficient scheme 
for implementation of two-qubit nonconventional geometric quantum gates, based on a dissipative large-detuning 
interaction of two-three-level atoms with a cavity mode is initially in the vacuum state . Many generalizations have 
been proposed to the original definition [38-42]. The basic theory behind geometric phases and their importance in 
quantum theory have been reported in [43l - l46j . In (47| experiments are proposed for the observation of the nonlinear- 
ity of the Pancharatnam phase with a Michelson interferometer. Jordan [48| has used the ideas of Pancharatnam to 
obtain a definition of phase change for partial cycles. The ideas of Pancharatnam were also used [H, [5(| to show that 
for the appearance of Berrys phase the system needs to be neither unitary nor cyclic [H, HI] , and may be interrupted 
by quantum measurements. However, this is not every thing, the topics on geometric phases and time-dependent 
quantum syste ms 1511 still attract extensive attention of a large number of investigators in various fields, including 
quantum optics |53l.l54{. condensed matter physics [30|,[55j], nuclear physics [56[, gravity theory [57j as well as molecular 
physics (molecular chemical reaction) (58|. To the best of our knowledge, the study of the topological Berry phase 
[21| of an atom-field interaction, within the framework of exact quantum electrodynamics in dispersing and absorbing 
media, has so far never been considered in references. A much richer range of phenomena is to be expected when 
allowing for the presence of dispersing and absorbing media, where a complex interplay of the electric properties of the 
atom and the bodies influences the functional dependence of the atom-field interaction. This problem is addressed in 
the current work, where we shed light on the topological Berry phase within the framework of the cavity QED, which 
is an important solid-state system for implementing quantum computation. We expect that, the additional noise due 
to the presence of dielectric media will, considerably, affect the measured values and behavior of the argument overlap 
between the wavefunction evolution and its original state. This approach, which has the advantage of being simple 
and applicable to different configurations of three-level systems, renders general expressions for the three-atom wave 
vector, as shown in Sees. II and III. In Sec. IV, application to the general results are examined for an atom placed 
in free space as well as in front of a dielectric half-space. A definition to the topological Berry phase and the form we 
are going to use including our numerical calculation and discussion are given in Sec. V. A summary is given in Sec. 
VI. 



II. GENERAL FORMALISM 



Consider a neutral atomic system with its center-of-mass positioned at the point ya near an infinitely long single- 
wall of dielectric surface. The total nonrelativistic minimal-coupling Hamiltonian of an atomic system interacting 



with a medium-assisted vacuum electromagnetic field in the form [61 



H = J d 3 rj^ doj hw f' f (r,u;) f'(r,w) 
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where rnj, Qj, Tj and Pj are, respectively, the mass, charge, position operator (relative to r^) and momentum of 
the jth particle constituting the atomic subsystem. In the Hamiltonian (pi, the electric field strength is expressed 
in terms of a continuum set of bosonic fields f(r,w) and f^(r, oj) which play the role of the fundamental (dynamical) 
variables of the composed system (electromagnetic field and the medium including the dissipative system) and satisfy 
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the well-known commutation relations 

[f m (r, W ),ft(r>')] =S mn 5{uj-Lo')5{v-v'), (2) 

[f m (r,o;),f n (r',a;0] = [ftfow),^.*/)] = 0. (3) 
In terms of the variables f(r, lu) and P(r,u) the medium-assisted electric field operator E(r) can be expressed as 

E(r)=P(r) + EH( r ), (4) 

with 

/>oo 

E (+) (r)=/ duE(r,w), E ( - } (r) = [E (+) (r)] t , (5) 

Jo 



E(r,w)=i,/A^ f dVv^MGtr.rVj.ftrV), (6) 

V £07T C Z J 

where G(r,r',w) is the classical Green tensor satisfying the equation 

TUJ 2 1 

-5-e(r,w) - V x V x G(r,r',uj) = -S(r - r'), (7) 

and satisfies the boundary condition at infinity, i. e., 

G(r,r',u) ->■ if |r-r'|->oo, (8) 

with complex permittivity e(r, w) is a function of frequency and space with the real part er{t,u)) and the imaginary 
part £j(t,uj) satisfy the Kramers-Kronig relations (for any r). The second and third terms of the Hamiltonian ((T|), 
represent the kinetic energy of the charged particles and their mutual Coulomb interaction, respectively, 

is the scalar potential of the charged particles distributed with the density, 

^(r)=^Q^(r-r i ), (10) 

3 

in the atomic subsystem. The last term of the Hamiltonian accounts for the Coulomb interaction of the particles with 
the the medium. 

A. Hamiltonian in the electric dipole approximations 

Applying the electric-dipole approximation, we may write the Hamiltonian (JT|) in the form (see Appendix A) 

//■oo 
d 3 rj du huP{r,w)-{(r,u)+Y^hu k S A kk-d A --E{r A ). (11) 

In this equation, the second term is the Hamiltonian of fc-level atom, where the SUfcfc are the atomic flip operators, 
while the last term is the atom-field coupling energy, where d A is the electric dipole atomic operator, [see Appendix 
A, Eqs. (|A9[ IA1 . where the medium-assisted electric field E(r) is expressed in terms of the variables f'(r,w) and 
f(r,w) according to Eas.(|5KH|). 
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FIG. 1: (a) Scheme of relevant energy levels and transitions for the three-level A-type atom, (b) An atom interacting with a 
half-space dielectric plate with a chosen atomic position za- 

B. Dynamics of a pumped three-level A-type atom 

Let us now focus on a three-level A-type atom, [Fig. 1(a)], near a one-dimensional infinitely long single- wall dielectric 
layer - the atom is on the right of the plate, with za being the (positive) distance between the surface of the plate 
and the atom [see Fig. 1(b)]. In particular, we assume that the atomic transition | 2) — >\ 3) is strongly coupled to 
the field modes via the dipole CI23. Moreover, we assume that an external (classical) pump field, with frequency uj pnl p 
and intensity described by the Rabi frequency Q, pmp , is applied to the atomic transition | 1) — >| 2). We may start 
from the one-dimensional version of the Hamiltonian (|11[) and apply the rotating- wave approximation [see Appendix 
A], the Hamiltonian (ITT]) reads 

H = J dz J™ cLj huf^r,u)f(r,u) + hu S 22 -(d 23 S 23 ^ (12) 

III. EQUATION OF MOTION 

In what follows we assume that the atom is initially (at time t = 0) prepared in a superposition of its two classically 
pumped levels | 1) and | 2), and the rest of the system, that consists of the electromagnetic field and the cavity, is in 
ground (vacuum) state | {0}). For simplicity, and without loss of generality, we assume that the frequencies of the 
atomic transitions | 2) o| 1) and | 2) 3) are the same; oj 2 i = ^23 = coo- We may therfore expand the state vector 
of the overall system at later time t > as 

//■oo 
dz J du C 3 (z,u,t) e-^ I 3) ft( ZjW ) I {o}), (13) 

where tf(z,w) | {0}) is a single-quantum excited state of the combined field-cavity system. It is easy to prove that 
the Schrodinger equation for | ip(t)) leads to the following system of (integro-)diffcrential equations for the probability 
amplitudes Ci(t), C 2 {t) and 6*3(2, uj,t): 

C l(t ) = j!W e iV P e lA ^C 2 (t), (14) 

with 
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C 3 (z,u,t) = -^tL^i^S) G*(z A ,z,u;) C 2 (t) e^-^K (17) 
Substituting the formal solution of Eq. (IT7|) [with the initial condition C^(z, 10, 0) = 0] 

C a (z,u,t) = -%!LJ^y/e>'{z,w) G*(z a ,z,uj) f dt' C 2 (t') e^-*)*', (18) 



\/Heo7r c : 

and of Eq. (fl"4"|) [with the initial condition C\{t = 0) = Ci(0)] into Eq.(|16p. and employing the integral relation 

uj 2 f 

— dze"(z,uj) G(z a ,z,lo)G*{z a ,z,lo) = Im G(z A , z A ,u), (19) 

we arrive at 

C 2 (t) = . ^^V '^ e-^'dfJ) + f dt' K(t,t') C 2 (t'), (20) 
2 Jo 

with the kernel function at the position z A 

\d I 2 P°° 

K(t,t') = -.p^ du uj 2 lmG{z A ,z A ,uj) e-^-"^-^. (21) 

It is worth noting that all the matter parameters that are relevant for the atomic evolution are contained, via the 
Green tensor, in the kernel function (1211) . 

IV. APPLICATIONS 

A. Three-level atom in free space near a perfectly reflecting mirror 

Besides being an interesting case in its own right, the free-space case will be useful in the comparison and interpre- 
tation of the outcomes. Since the electromagnetic vacuum cannot be switched off, its interaction with atomic systems 
cannot be switched off either, thereby giving rise to a number of observable effects such as spontaneous emission, the 
Lamb shift, intermolecular energy transfer, and so forth. In free space, for a sufficiently high-Q vacuum cavity con- 
tains nondispersing and nonabsorbing media, the excitation spectrum of the electromagnetic field effectively extends 
over a frequency interval, so that the w integrals [see, e.g., (|21[) ] effectively run only over this interval. If we assume 
that the frequency interval is sufficiently small and sufficiently far from medium resonances, so that both dispersion 
and absorption may be disregarded. In this case, the permittivity becomes approximately real and independent of 
frequency. Note that a constant and real permittivity in the whole frequency domain can be realized only for unity 
permittivity, i.e, the vacuum. This can be used to express the vector potential and the canonical momentum field 
in terms of the monochromatic modes and the associated photon annihilation and creation operators. In short, the 
excitation spectrum turns into a quasi-discrete set of lines of mid- frequencies u>k and widths Tk- Furthermore, in 
the emission spectrum, the intensity of light associated with the emitted radiation is going to be centered about the 
atomic transition frequency loo, with widths Tk ~ 0. Bearing in mind these physical interpretations, with the Green 
tensor G(z A , z a ,uj) is replaced by the vacuum Green tensor; G v (z A , z A , u>q) [61J, where 



JmG v (z A ,ZA,uo) = ^l, (22) 

D7TC 



equation (|2"Tj) becomes 



K (t,H) = -^Lp- I du k col e'^-^-t'l (23) 
neoirc^ one J 



The quantity uif, varies little around uj^ ~ ujq, we can therefore replace uSf, by ojq. The integral 

duj k e -i(«*-"o)(t-t') = 2n5(t - t'). (24) 
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yields the following equation for C2(i), in the Weisskopf-Wigner approximation 



C 2 (t) = %. 



<*pmp -i$„ 



-iA„t 



»*Ci(t).-r^c 2 (t) 



with 



fel^ 3 



p32 _ 

' 3fen7TC 3 



is the free-space spontaneous emission rate for the transition | 2) — >| 3). 

Upon taking the Laplace transforms of equations (|T4"|) and (|25|) with the use of the initial conditions [C\ (t 
Ci(0) and C 2 (t = 0) = C 2 (0)], we find that 



C 3 (*) 



(a + zA p )C 2 (0) + jfa e -^ Cl ( ) 
p( s ) 



and 



Ci(s + iA p ) = 

where Cj are the Laplace transforms of Cj, and 

p{s) = s 2 + 



[s + rg 2 ]Ci(0) + z%^e"V^C 2 (0) 



rf + iA p 



fe=i 



and Xk, (k = 1, 2) are the roots of the quadratic equation 



rf + zA p 



ft 2 

V"P _i_ ,t 32 A 



found by substituting x = s in the equation p(s) = 0, they are given by 



r 32 + zA r 



£1,2 



± 



rr 32 - zA p 



2 n 2 

pmp 



If the classical pump laser field is resonance with the atomic transition frequency; A p — 0, we have 



-:S2 



£1,2 



± 



(r 32 ) 2 n 2 p 



pmp 



2 V 4 4 
Using the residue theorem, Eqs (|27|) and (f28| are easily inverted to give 



Ca(*) = E 



8<t (jEj + iA p )(7 2 (0)+i 



g LfJ pmp 



d(0) 



(25) 

(26) 
0) = 
(27) 

(28) 

(29) 

(30) 
(31) 
(32) 
(33) 



e-^dW^Ee- ^'^^^'T^^ 2 ^ U = M; Mj (34) 

with 

d(0) = -^e* C a (0) = Ci(0)e-^ (35) 

From equation (1321) we sec that both roots lEj are (a) negative when fl pmp ^ Tq 2 , and (b) complex (with a negative 
real part equal to — Fq 2 /2) when O pmp > Tq 2 . Thus the time evolution of amplitudes C m (t),m = 1,2 (and hence of 
the upper- level populations can be divided into two regimes of different behavior. For tt pm p > Tq 2 , the populations 
display pronounced oscillations before decaying to zero. On the other hand, when £l pmp < Tq 2 the populations barely 
complete an oscillation before decaying to zero [62| . Thus, the driving field induces oscillations on the populations of 
the upper levels. The stronger the driving field (i.e. the larger the ), the faster the oscillations. 
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1 . Driving field with intensity equals to the vacuum spontaneous emission rate 

Thus, as we seek for comparison, we consider only the case when the pump laser field is applied with intensity; 
Qpmp, equals to the free space spontaneous emission rate Tq 2 , i.e., 0, pmp = Tq 2 . When this occurs, equation (|3"2"j) has 
a double root Xi = x 2 = —Tq 2 /2 and inversion of equations (|3"51) and gives 

r 32 r V 32 t / \1 

C 2 (t) = e— 8-*[c a (0) - -|^(c 2 (0) - «>-**"»* Ci(0))] (36) 
r 32 r r^ 2 ^ / \ i 

C x {t) = e"^ 4 [d(0) + -|-(Ci(0) + ie**™ C 2 (0))] , (37) 
which in agreement with results of Ref. [62] . 

B. Three-level atom in front of a perfectly conducting planar dielectric half-space 

The presence of linear media in the form of macroscopic bodies changes the structure of the electromagnetic field 
compared to that in the free space and in consequence the electromagnetic vacuum felt by an atom is changed. As we 
mentioned above, absorption is always associated with additional noise which its effect included in the corresponding 
dynamical variables of the medium. In addition, the photonic density of states can be modified by the presence of 
macroscopic bodies, accordingly the mathematical treatment of the integral kernel given by Eq. (121 [) will be completely 
modified, and as a consequence its physical interpretation. Starting from Eq. (|20[) . on taking the time integral of 
both sides of this equation, it is not difficult to see that the equation can be converted into the integral equation 

C 2 (t) = C 2 (0) + !fe e -»-p /' dt' e-^'dit') + f alt' K{t,t')C 2 {t') (38) 
^ Jo Jo 



K(t,Q = r^^ Im f*^"> ( e-^Kt-O _ i ) (39) 



where 

k ■ i /' , 

Heott J c 2 i(uj — ujq) 

From the above it is seen that C 2 {t), with the presence of macroscopic bodies being fully included in the Green tensor 
of the system, depends crucially on the actual structure of the Green tensor which carry all information about the 
macroscopic bodies properties. In general, it seems not easy to calculate C 2 (t) in an exact analytical form. However, 
using limiting cases of Eq. (|39D will simplify the problem to obtain the analytical solution without explicitly making 
use of the actual structure of the Green tensor. 

In order to study the case where the atom is surrounded by matter, the atom should be assumed to be localized in a 
small free region, so that the Green tensor at the position of the atom reads 

G{z A , z a ,uj) = G v '(za,Za,u) + G r (z Ai z a ,uj) (40) 

where G R (z Al z A , ui) describes the effect of reflection at the (surface of discontiniouty of the surrounding medium), 
and G v (z A , z A ,u) is the vacuum Green tensor given by Eq. (l22l) . The contribution of G v to K,(t, t'), in Eq. ([39]) takes 
the form 

jc(t,0 = -!*»!! r d ^ G J^ z ^( e ^'^ l) 

He iT J c l i(u - u) \ J 

l^ 23 ' 2 f°°duj- Im GR ^ Z ^) ( c i(u -u)(t-f) A (41) 



This equation can be regarded as the basic equation for studying the influence of an arbitrary configuration of 
dispersing and absorbing matter on the spontaneous emission decay and related phenomena. 
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B.l. Weak atom-field coupling 



To achieve our goal, and hence, demonstrate the usefulness of this article simply, let us restrict our attention 
to the case when atom-field system is coupled weakly. This can be done mathematically by applying the Markov 
approximation. When this approximation applies, i.e., when in a coarse-grained description of the atomic motion 
memory effects are disregarded, then we may let 



i{wa— w)(t— t') _ i 
i(ui - uj) 



in Eq.gU) [((x) = %S(x) + iV\], and thus 



fC{t,t') 



p32 
"~2~ 



iSujQ 



where T 32 and Suiq are respectively given by 



(42) 



(43) 



(44) 



SuJo = l^L V l ^ ^G(z a ,z a ,uj) 



(u - UJ ) 



(45) 



Obviousely, in contrast to free space, the decay rate now becomes a function of the atomic position and the orientation 
of the transition dipole moment. In other words, the vacuum fluctuations felt by an atom in the strictly free space are 
inhomogeneously and anisotropically changed by the presence of the bodies, in general. The nearer to a body that 
the atom is located, the stronger the effect to be expected. 

Following the approach of Dung et al., [f|, and after Recalling the Kramers-Kronig relation for the Green tensor, 
we may approximately rewrite Suiq as 



Sujq 



Hsqc 2 



Re G (z A , z A ,u a ) 



duj 



uj 2 Im G R (z A , z A , uj) 
ul (uj + w ) 



(46) 



The second term, which is only weakly sensitive to the atomic transition frequency, is small compared to the first one 
and can therefore be generally neglected. Obviously, the divergent contribution of vacuum to the line shift is included 
in the atomic transition frequency ujq. Subistituting (|43[) into (|38[) and taking the time derivative of both sides with 
the setting of T 32 /2 = ot 32 for simplicity, we obtain 



c 2 (t) = jM e -i<„ e -iA,./ 



Ci(t) + (-vo 32 + iSujo) C 2 (t), 



(47) 



the next step is to apply the Laplace transforms on Eqs. (IT41 and (|4"7]). with the same initial conditions, to obtain 



and 



where 



C 2 (s) = 



d(s + iAp) 



(s + iA p )C 2 (0) + l hl^ i p e -^ PmP Cl (0) 

pO) 



s - (-to 32 + iSuo]Ci(0) + V-^ e ^^C 2 {Q) 



(48) 



(49) 



s) = s z + [pj iz - i5uj ) + iA p s + iA p (vj iz - iSuj ) + 



pin p 



k=l 



(50) 



with yj~ are given by 



2/1,2 



(tn 32 — iSoJo) + iA p 



± 



(tu 32 — iSuio) — iAp 



2 Q 2 

pmp 



(51) 
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Following the same procedure, as previous, it is not difficult to obtain C m (t), m = 1, 2 in the form 



(2/i - 2/j) 
- (-tn 32 + i5wo)]Ci(0) + 



B.l.l. Model permittivity of Drude-Lorentz type 



e -^C 2 (0) . . 



(52) 



(53) 



To be more specific, let the classically-pumped three-level atom located near the surface of an infinite half-space 
dielectric medium such that 



e(r,w) 



if z < 
1 if 2 > 



(54) 



To give an impression of what can be observed in real situation, let choose a planar dielectric surface of permittivity 
modeled by the widely-used-in-practice Lorentz type as 



e(w) = 1 + 



iW7 



(55) 



where lot and 7 are the medium oscillation frequencies and linewidths, respectively, and lo p correspond to the coupling 
constants. An example of the permittivity for a (single resonance) dielectric as a function of frequency is shown in 
Fig. 3. For z > 0, and when the distance of the atom from the surface is small compared with the wave length, 



£/ J 









0.90 0.95 



1.15 




1.06 1.08 



1.00 1.05 1.10 

FIG. 2: (a) Real and imaginary parts of the permittivity of single-resonance Drude-Lorentz type dielectric for uj p = 0.5ujt, 
7 = 10 _4 cjt- The band gap covers the interval from ujl = Ut to uil — 1.12ujt, (b) Spontaneous decay rate T 32 [Eq. <|56|) ] and 
level shift 5uo [Eq. (|57|) ] for same lj p while 7 = 10~ 2 cjt and za = 0.05At. 



fcz< 1, using the results obtained by Scheel et al. [60( of the reflection part of the Green tensor yields for the decay 
rate 



^32 



3r 3 



:S2 



c 

LV Z 



£j(^o) 
kM + l| s 



(56) 



The leading term (1=3 z~ 3 ), which is proportional to et(u}q), is closely related to the virtual photon emission with 
subsequent medium quasiparticle excitation (nonradiative decay), i.e., energy transfer from the atom to the medium. 
Similarly, for the line shift, due to the presence of the macroscopic bodies reads 



5u) 



32 



1 + — 

d? 



c 

UJQZ 



kM + i| s 



(57) 
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V. THE GEOMETRIC PHASE 



If the final wavefunction cannot be obtained from the initial wavefunction by a multiplication with a complex 
number, the initial and final states are distinct and the evolution is noncyclic. The total phase has both dynamical, 
4>d, and geometric phase, 4>g-> parts for an arbitrary quantum evolution of the system from an original state at t = 
to a final wavefunction at time t. Suppose state | ip(0)) evolves to a state | ip(t)) after a certain time t. If the scalar 
product [3!| 

(V(0) \exp^J H(l/)dt') | V(0)> (58) 

can be written as j3 e" 7 , where j3 is a real number, then the noncyclic phase due to the evolution from | ip(Q)) to 
| tp(t)) is the angle r\. This noncyclic phase generalizes the cyclic geometric phase since the latter can be regarded as a 
special case of the former in which (3 = 1. The total phase (the topological Berry phase) acquired during an arbitrary 
evolution of a wavefunction from the vector | ip(0)) to | tp(t)) is prescribed by Phancharatnam as 

fa = arg(VXO) | ^(0)> (59) 

On substraction the dynamical phase 4>d from the total phase 4>t we obtain the geometric phase <p g . An exact expression 
of the geometric phase can be obtained as 

& = arg{CJ(0)Ci(t) + C 2 *(0)C 2 (i)e- lWo *} 



arcsin ( Y ^ ) , (60) 



\y/X*(t)+Y*{t). 

with 

X{t) = Re[Ct{0)Ct(t) +C* 2 (0)C 2 {t)e- Wot ] 1 (61) 

and 

Y{t) = Im(Cr(0)Ci(t) + C 2 *(0)C 2 (i)e-^ ot ), (62) 

where the dynamical phase for an arbitrary quantum evolution from a time t — to a time t is given by the time 
integral of the expectation value of the Hamiltonian over the time interval [0, t] 

m) | H{t) | ij{t))dt (63) 

As might be expected, the behavior of the present three-level system changes dramatically depending on the variation 
of the adjustable parameters of the system, namely, dipole vector orientation, the atom position, za from the surface, 
and finally, the surface permittivity. 




FIG. 3: Two dimension plot of | 2} state population (a), the topological Berry phase <f>t (b) and the fluctuations in X and Y 
for atom in free space with Qprnp/To 2 — 1. 



Examples [for all our plots we have used fixed values for Ci(0) = (^(O) = -75 and cj) — n/2] of the populations and the 
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corresponding topological Berry phase against Tg 2 i are shown in Fig. 3 for an atom in free space. It is rather interesting 
to mention to the fact that: in the Feynman approach to quantum mechanics it is important not only the modulus 
and argument of complex number, but the quantum fluctuations of this variables too. In such a representation, 
the deviation of quantum trajectories from the classical trajectories are described by the quantum fluctuations (63[ 
. Therefore we have further investigated the fluctuations of the real and imaginary parts of a.rg(ip(t) | "4>{t)). At 
first, for initial C2(0) = 0, the topological Berry phase <pt vanishes, this means physically that there is no phase in 
photon transition that make the photon transition is in a straight line, i.e., the photon transition will occurs without 
any phase. Thus the atomic coherence is necessary for the development of the topological Berry phase <j) t . Note the 
populations decay exponentially while <pt normally oscillates, due to the oscillations in X and Y, with an amplitude 
increases with time. 

In an interesting proposal for detecting the topological Berrys phase was suggested [64], (65[ for a certain adiabatic 
evolution of a joint state of the internal levels of a trapped ion and its vibrational motion. It has been shown that 
despite being the phase gained by a joint state, its value is fundamentally dependent on the harmonic oscillator nature 
of the vibrational mode. 

When the effect of dielectric matter is considered, we clearly notice very different figure shape. It worth to note 
that, in the short range interval of ujq/uit at which the decay rate makes half period and change its evolution direction, 
the level shift 5ujq sa 0, in this interval the complex amplitudes C m (t), m = 1,2, Eqs. (|52[ [53l) . after putting A p = 0, 
depend only on the decay rate T 32 , and can take the same form of Eqs. (|3"6l [37]) of free space with the replacement of 
the free space decay rate Tq 2 by w 32 . Examples of the population, (f> t and the fluctuations in X and Y for dielectric 




FIG. 4: Mesh plot of | 2) state population (a) and the topological Berry phase <j> t (b) as functions of the transition frequency 
loo/lot and F^t near a planar dielectric half-space for an z-oriented transition dipole moment, with f2 pmp = to 32 , lo p = 0.5ujt, 
z = 0.05A T and 7 = 10~ 3 lj t . 

matter of Drude-Lorentz type, with linewidth 7 = 10~ 3 , are shown as mesh plot in Figs. 4, 6 and 8 with the 
corresponding cross section (the dependance on the transition frequency) when Tfft — 10.0 in Figs. 5, 7 and 9. The 
property of Drude-Lorentze type, where it features band gab between the transverse frequency lot and longitudinal 

frequency lot = J w T + w 2 , where Re e(w) < 0, control strongly the behavior of the functions under study. Moreover, 
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FIG. 5: Cross section of Fig. 4 as well as fluctuations in X and Y (c), when Iq t = 10.0 while rest of parameters as in Fig. 4. 



since the Drude-Lorentze model incorporates surface-guided waves, which are observed in the presence of the interface 
and bound to it, with the amplitudes damped into either of the neighboring media, we notice the populations, <p t 
and the fluctuations in X and Y affected strongly with these properties and appear as surface-guided functions. It 
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is seen that, in the band-gap region, due to the drastically increase if of the spontaneous decay rate, because of the 
non-radiative decay channel associated with absorption, the population amplitude damped into the band-gap region 
while (f> t approaches unity as seen from the figures. Fig. 5 also shows that the populations vanishes inside the band 
gap in the interval where ujl S [wti 1.12ojt], as expected, while fluctuations go beyond negative values. On increasing 




FIG. 6: The same as Fig. 4 but for 7 = 10~ 2 Wt 

the damping parameter^, one observes, the narrow-band gradually getting so much broaden [Figs. 7; 7 = 1CP 2 and 
9; 7 = 10 _1 ] that it eventually disappears (not shown), due to the increase of spontaneous decay for frequency ojo 
- with value equals to unity for <f> t is preserved- where most of the energy emitted by the atom is absorbed by the 
medium wall in the course of time while a fraction of light can scape to free space which reflects the clearly-noticed 
decay of the population as time passes [see Figs. 6(a) and 8(a)]. Moreover, in the absorption band and in the vicinity 

of the absorption band, i.e., in the region between resonance uit and longitudinal frequency uil = *J ^\ + ^>p (in the 

figures, ujl = 1.12WT), these functions respond sensitively to the change of the damping parameter 7. The differences 
between the cases here are less pronounced for strong absorption, i.e., when the value of the bandwidth parameter 7 
is sufficiently large (compare Figs. 5, 7 and 9). We conclude that the differences obtained differ essentially in the way 
the bandwidth parameter is introduced. The maxima of the populations below the band gap and at the field resonance 




(cjo ~ wt) frequencies (Figs. 4 and 5) indicate that, although the decay can be noticeably enhanced, the probability 
of emission of a really observable photon can be substantially reduced compared to the case of spontaneous emission 
in the free space. Obviously, a photon emitted at such a frequency is typically captured by the surface for some time 
due to the deep penetration of radiation into the surface, and hence the probability of photon absorption is increased. 
Also, for transition frequencies inside the band gap, two regions can be distinguished. In the low-frequency region 
(cjo < ut), where surface-guided waves are typically excited, radiative decay dominates, i. e., the atomic transition is 
accompanied by the emission of a photon escaping from the system. Here, the radiation penetration depth into the 
surface is small and the probability of a photon being absorbed is also small. With increasing of atomic transition 
frequency the penetration depth increases and the chance of a photon to escape drastically diminishes. As a result, 
photon absorption dominates [68j . In this case, in the regions below and above the band gap, 4>t preserves a fixed 
negative short angle with |</>t| ~ 77°, which means that the wavefunction evolution and the original state move apart 
from each other, while inside the band gap region this angle decreases to reach a short angle equals to <fi t ~ 52° 
(see Fig. 5). On increasing 7, 4> t differs considerably in the two regions, where below band gap <p t shiftes to become 
smaller than above band gap region (Figs. 5 and 9) due to the strong penetration of the emitted photon into the 
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surface. We can say that the topological Berry phase strongly depends on the presence of dielectric surface and such 
systems are potentially interesting for their ability to process information in a novel way and might hnd application in 
models of quantum logic gates. Figures 10-15 illustrate the dependence of the populations, topological Berry phase, 




FIG. 10: Fhe same as Fig. 4 but when taking into account the effect of both T as well as Scoo 

4>t > and the fluctuations in X and Y on line shift 5u>o as well as the decay rate Tq 2 of an atom surrounded by a planar 
dielectric half-space with various values of the bandwidth parameter 7. The figures reveal that with the decrease of 
bandwidth parameter 7, specially outside the band gap region, where ; loq < ujt, i.e., s(u>o) > 0, the surface-guided 
waves, due to Drude-Lorentz model, give rise to strong enhancement of functions' oscillations to the extent that the 
lines strongly overlap in all time stages regardless of the population decay. The same effect can be easily seen in X 
and Y with particular emphasis on that X reveals oscillations in all positive half plane outside the band gap region 
while goes once to the negative half plane inside the band gap region before accomplishes its decay. The function Y 
oscillates around zero with approximately the same behavior of X. Also, for sufficiently small values of the bandwidth 
parameter, 7, below the band gap, energy transfer between the atom and the medium occurs rapidly to the extent 
that the populations of the levels | 2) and | 1) overlap except for small course of time, while energy transfer strongly 
enhanced above the band gap region, which clarifies the high maxima and minima of (fit [see Fig. 11]. On increasing 
7, populations overlap breaks down gradually and, at end, energy transfer stops with high values of 7, where the 
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atom preserves its lower state and energy transfer occurs only above band gap region with small rate as can be seen 
by comparing Figs. 13 and 15. This means a small separation angles can be produced. All these results support our 
idea that selection of the atomic state leads to a higher degree of nonclassical effects in the resulting radiation field. 
Now we turn our attention to the effect of the distance za of the atom from the surface of the plate. Examples 





1.0 1.5 



1.0 

oj /a> T 



FIG. 13: The same as Fig. 11 but for 7 



of the populations, topological Berry phase, (f>t, and fluctuations in X, and Y, as functions of za in the vicinity of 
a surface-guided field are shown in Figs. 16 - 21, where a single-resonance Drude-Lorentz-type dielectric has been 
assumed. One can see that the oscillations collapse after a few Rabi periods - for small interval below the band gap- 
and after an interval of time topological Berry phase 4>t vanishes, which means an in phase photon transition. In this 
case, an energy transfer between the atom and the surface occurs for distances short from the surface leaving the atom 
dexcited regardless of its distance from the surface due to the deep penetration of the emitted photon into the surface. 
This means, physically, that, in this case, the atomic transition frequency is relatively near medium resonance, so that 
the imaginary part of the permittivity becomes relatively big which enhances medium absorption. In addition, inside 
the band gap region, (j) t preserves fixed value - even when energy transfer occurs for sufficiently small values of 7 [Fig. 
16] - regardless of the atom being far or close to the surface plate. As 7 increases Rabi periods shortened gradually 
where the general behavior preserved. So one can clearly see that the topological phase 4> t is affected strongly by 
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the presence of the half space for atom-half-space separations that are greater than the medium oscillation wave 
length At (</>t ~ for za ^ 0.5At as the figures suggest). These results means that by adjusting the bandwidth 
parameter 7, Rabi strength of the driving classical, ^l pmp and the distance between the atom and the surface, Za, 
quantum mechanically, we can with full success, through adiabatic quantum process, expect that the cyclic evolution 
of wavefunction yields the original state plus a zero phase shift, namely, cyclic evolution of wavefunction develops on 
the same track of the original state. In other words, the angular separation between the wavefunction evolution and 
the original state dies, that is of strong applications in construction of the universal quantum logic gates. 

In fact, in quantum computation, operations are performed by means of single-qubit and multiple-qubit quantum 
logic gates [H, [fill- One could use the present model to generate a C-NOT gate, i.e., universal quantum logic gate 
based on three-level atom. For polarising beam splitter (PBS), for an emitted photons, PBS redirects vertically 
polarized photons ( in state | V)) without affecting horizontally polarized photons (in state | H)) such that a C-NOT 
gate can be obtained as 

a\H) 1 +/3\V) 2 ^a\H) v +(3\V) 2 ,, (64) 
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FIG. 17: Cross section of Fig. 16 as well as fluctuations in X and Y (c), when Vq t = 2.0, tv — 0.5lj t and 7 = 10 3 iv T , 



where, for phase plate (PP) which adds an overall phase factor to the state of the incoming photon such that 

I h. e** I (65) 

with phase factor <fi = 0, for present case. This can be accomplished my moving the atom away from the surface 
of the half-space. Note that there are values of the atom-surface separation at which the topological Berry phase is 
typically trivial [Figs. 17 and 18]. 




FIG. 19: The same as Fig. 17 but for 7 = 10" 2 o;t. 



VI. SUMMARY AND CONCLUSIONS 

Within the frame of macroscopic QED in dispersing and absorbing media, we have investigated the topological 
Berry phase properties of the photon transition between two classically pumped levels of three-level atom interacting 
with the electromagnetic vacuum field in the presence of dispersing and absorbing material surroundings described by 
a spatially varying permittivity that is a complex function of frequency Using a source-quantity representation of the 
electromagnetic field, and starting from the exact minimal-coupling Hamiltonian (in the electric dipole approximation) , 
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FIG. 21: The same as Fig. 17 but for 7 = 10" 



we have performed the calculations for a one-dimensional Hamiltonian of the total system that is accordingly used 
in the calculation of the wavefunction, in a general form, that governs the system utilizing Schrodinger equation. To 
gain further insight into the difference in the nature of the topological Berry phase function, we presented a more 
subtle comparison to include the atom in free space. Without ab initio calculations, utilizing the physical properties of 
vacuum field evolution applied to the kernel function K (t, t'), we have derived the complex amplitudes of the classically 
pumped levels in free space, hence, the argument overlap between photon transition is also calculated applying the 
Weisskopf-Wigner approximation. On considering the atom near a dielectric body, we restricted our selves to a 
planar dielectric half-space with model permittivity of Lorentz type. Effects of various physical settings regarding the 
transition frequency, bandwidth parameter and the position of the atom, za, relative to the dielectric surface have 
been investigated. The results show that, in contrast to free space case in which the probability amplitudes of the 
classically pumped levels decay exponentially, the emitted light observed, and the topological Berry phase between 
photon transition can be controlled strongly by the presence of macroscopic bodies as follows: 

• For relatively fixed small distance za, between the atom from and dielectric surface, specially, in the case 
that the line shift is canceled, the topological Berry phase exhibits maximum angles separation between the 
wavefunction evolution and original state in two positions, namely, before interaction is switched on and in 
the band gap region where the decay rate reaches its maximum. This behavior dominates accompanied with 
band-gap-envelope phase gets more broaden as bandwidth grows. When line shift is taken into consideration, 
the phase shows Rabi oscillations, specially below the band gap region, while these oscillations invoked below 
and above band gap with bandwidth increases. 

• Moving the atom a way from the dielectric surface leads to phase angle evolution in straight line coincide with 
x-axis, which means an in-phase photon transition frequency. 

• For small separation distance for atom placed in parallel alignment with respect to a half space, the dielectric 
surface-induced modification of the topological Berry phase is very different from our findings of atom placed 
in free space: A dielectric surface reduces the angle separation between the wavefunction evolution and the 
original state, while free space leads to an enhancement. In other words, if the atom is placed sufficiently 
near the dielectric surface, its spontaneous decay may be suppressed, with the emission line being accordingly 
narrowed, which encourage energy transfer between the atom and the dielectric medium. 

• Further, in contrast to a classically pumped atom in free space, the exponential decay dynamics gives place to 
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the non-exponential one if the atom moves away from the dielectric surface due to the coexistence effects of the 
dielectric body and classical pumping laser field. Thus, the atom-field coupling strength and the character of 
the spontaneous decay dynamics, respectively, may be controlled by changing the distance between the atom 
and the dielectric surface by means of a proper preparation of classical pumping laser field. 

• At end, it should be stressed that transition frequency shifts experience significantly enhanced changes than in 
free space case. 

• Our study opens routes for new challenging applications of atomically systems as various sources of coherent 
light emitted by pumped atoms in dielectric surroundings. 

• The (numerical) calculations performed for atoms near a half-space can of course be extended to other geometries 
of the material surroundings such as microspere or cylindrical objects. 



Appendix A: THE HAMILTONIAN IN THE DIPOLE AND ROTATING WAVE APPROXIMATIONS 

Applying the minimal-coupling scheme (flj , we may decompose the Hamiltonian (jTJ) of the coupled system consisting 
of the atom and medium-assisted electromagnetic field as 

H = H F + H A + H AF (Al) 

where 

H F =Jd 3 rJ^ doj fiw f /f (r,a;) f'(r,o;) (A2) 
is the Hamiltonian of the electromagnetic field and the dielectric matter, 

= it- + \ I d3r P^A(r), (A3) 

is the Hamiltonian of the atom, and 

Haf = J2 —V ■ + / d 3 rp A (rMr), (A4) 

3 J 

is the interaction energy. Here we omitted the A 2 term as a result of the (quasi-)resonant atom-field interaction. In 
the electrc-dipole approximation the first term on the right-hand side of (|A4|) simplifies to 



- J2 — Pi ■ Mrs) = E ' Mr a) = -±[d A , Ha] ■ Mr a), (A5) 

Tit i % ih zrt 

j 3 

where 



d A = Qjfo ~ TA )=Y1 $3*31 ( A6 ) 



is the atomic dipole operator for a neutral atom with the nucleus being positioned at ya- Note, in (| A5|) the dipole 
approximation has been employed by replacing A(r 3 ) — > A(rx) and the commutation relation pj = — ^-g 1 [r,- , Ha] , 
bearing in mind that in the Coulomb gauge [pj, A] = 0, has been used. 

In the next step, the atomic Hamiltonian (for an isolated atom) satisfies the eigenvalue equation 

H A I k) A = hw k | k) A , (A7) 

where Hulk are the energy eigenvalues and | k) are the eigenvectors. The states | fc) span the Hilbert space of the 
atomic system. Since Ha is an observable, its basis vectors {| k)} form a complete orthonormal set as expressed by 
the closure (completeness) and orthonormality relations 

Y,\QAA(k\=I, A(k\k') A = 5 kk >, (A8) 
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where I represents the identity operator, and 5kk' is the Kronecker delta function (Skw — 1 f° r k = k' , and Skk' = 
for k k' ). We dehne the atomic operator §Akk' by 

S Akk > =\ k) AA (k' | . (A9) 

Using Eq. (|A8[) . Ha can be written as 

H A = ^2 t^kSAkk- ( Al 0) 

k 

Similarly, on using Eq. (|A8j) twice, we can represent the atomic dipole-moment as 

d A = 22 Akk Akk '> (AH) 
k,k' 

with 

<W =a (fc | d A | (A12) 
If we apply the same procedure on Eq. (|A5p . we arrive at 

- j^\dA,H A ] ■ A(r A ) = iy)uk>kSAkk> d A k'k • A(r^); Wfc&' = w,t - = —a;*/*,, (A13) 

fe.fe' 

This equation applies to a general fc-level atom interacting with a quantized electromagnetic field in the electric dipole 
approximation. In the following we will restrict our attention to a three-level system. A three-level atom (consisting 
of the ground level | 3) and excited levels | 1) and | 2)) can be in one of three distinct ( cascade; S, lambda; A and 
V) configurations. Thus, for each of the three configurations there is a direct dipole transitions is forbidden and two 
dipole-allowed transitions. Restricting our attetion to A-type, direct dipole transition is forbidden between levels | 1) 
and | 3). With the property of d A kk' where, for states of the same parity dAkk — 0, bearing in mind that dAkk' is 
Hermitian; dAkk 1 — d* Ak , k , equation (|A5j) . with use of Eq. (|A13j) . becomes 



SaiZ d,423 — SaZ2 dAZ2 



•A(r A ) (A14) 



In the Coulomb gauge, the vector and scaler potentials are given, respectively, by 

r°° 1 „ 

A(r) = / duj— E x (r,a;)+H.c (A15) 
Jo IUJ 

f°° 1 

- V<p(r) = EH (r) = / duj— E"(r,w) + H.c (A16) 
Jo lu 

where the transverse (longitudinal) electric field operator 

E- L <H)(r) = J dV^H>(r-r') •E(r'), (A17) 

with i5 _L (r) and <5"(r) being the transverse and longitudinal S- functions respectively. Recalling Eqs. (|A15I - ET7| . with 
the application of the rotating wave approximation, and after short calculations Eq. (|A5I) reads 



= -dA-E"(r A ) ~ -(s A23 M(+\r A ) ■ d A23 + H.c.) (A18) 

3 

where the total electric field is given by Eq. ((4]) with the transverse (longitudinal) electric field operator reads 

/>oo 

E x dD(r A )= / dwE-J-CID^.w) (A19) 
Jo 
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and ui = W23 is set in the integral. 

In order to deal with the second term in Hafi we expand Pa{ v ) in a multi-polar form and retain only the first 
non- vanishing term 



PA 



(r) = Y, QA* ~ ta) - V • (S(r - v A ) £ Qfa - r A 



= -V • S(r - r A )d A ^ - J d 3 r [v • 5{r - r A )d A ] <p(r) 

= J d 3 r[S{r - r A )d A ] ■ V0(r) = d A ■ V£{r A ) = -d A ■ E" (t a ) (A20) 

where integration by parts and Eq. (|A16[) have been employed for deriving the second and the third equation, 
respectively. Recalling Eq. (|A11[) . restricting our attention to A- type with the same properties applied to (IA14[) . we 
arrive at 

J d 3 rp A {r)£(r) = -d A ■ M(r A ) ~ -[S A23 M^ (v A ) ■ d A23 +K.C.} (A21) 

Combining (|A^| . (|A18|i . and (|A21[) gives 

H AF = -d A ■ E(r A ) ~ - [5a 23 E (+) (t a ) ■ d A23 + H.c] (A22) 
For three level atom, k = 1,2, 3, the atomic Hamiltonian (IA10[) reads 

H A = hu k S Ak k = fkjJxS A ii + hhj 2 S A2 2 + hw3S A33 (A23) 

Equations (|Alj) . (|A2|) . (|A22| . and a subtraction of ft(a;2 — u>o)I from ifA, Eq. (|A23j) (taking into consideration that 
CJ21 = ^23 = ^o), lead to the Hamiltonian 

H = J d 3 vj^ doj i'Hv.u) f'(r,w) + fiuoSjw - [S A23 ^+\r A ) ■ d A23 +R.c] (A24) 
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